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Abstract 

The bright soliton solutions of the mixed 2-coupled nonlinear Schrodinger (CNLS) equations 
with linear self and cross coupling terms have been obtained by identifying a transformation that 
transforms the corresponding equation to the integrable mixed 2-CNLS equations. The study on 
the collision dynamics of bright solitons shows that there exists periodic energy switching, due to 
the coupling terms. This periodic energy switching can be controlled by the new type of shape 
changing collisions of bright solitons arising in mixed 2-CNLS system, characterized by intensity 
redistribution, amplitude dependent phase shift and relative separation distance. We also point 
out that this system exhibits large periodic intensity switching even with very small linear self 
coupling strengths. 
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I. INTRODUCTION 



The study on soliton collisions has been receiving sustained attention since the identi- 
fication of their particle like collision behaviour by Zabusky and Kruskal [l| in 1965. Due 
to their intriguing collision properties and their robustness against external perturbations 
solitons find applications in diverse areas of science which encompass the current thrust 
research areas including nonlinear optics, Bose-Einstein condensates, plasma physics, bio- 
physics . Particularly, it has been shown that soliton propagation in systems such 

, multimode fibers |2], fiber couplers photorefractive medium 
[9|], Bose-Einstein condensates [lOj and continuum limit of Hub- 
bard model in ID is governed by multicomponent nonlinear Schro ding er type equations 
which become integrable for specific choices of system parameters 



as birefringent fibers [6 
left handed materials 



12 



urn- 



Recent theoretical and experimental studies show that the integrable CNLS system with 
focusing type nonlinearity exhibits fascinating shape-changing (energy /intensity redistribut- 
ing) collision of bright solitons characterized by intensity redistribution among the colliding 
solitons in all the components, as well as amplitude dependent phase-shifts and change in 
relative separations distances [15|, llg, |l7| • In such a two soliton collision process there occurs 



suppression/enhancement of intensity in few components and enhancement /suppression of 
intensity in the remaining components after interaction. We call such a collision scenario 
as type-I shape changing collisions (type-I SCC). This interesting collision behaviour has 
also been experimentally verified in birefringent fibers [3] and in photorefractive media 
[3]. Further investigations on this type-I SCC has led to exciting applications in collision 



based optical computing 



13, 



20. 



nonlinearities are of mixed type 



21 



23 



22l | . However the collision scenario is different if the 
24| , which includes both focusing and defocusing types. 
The corresponding mixed CNLS equations admit shape changing collision of bright solitons 
in a quite different manner from the collision scenario of type-I SCC. Very recently, it has 



been shown that in mixed CNLS equations during a two soliton collision process there is 
a possibinty of either enhancement or suppression of intensity in a given sohton in a., the 
components [24|. Here also the collision process is characterized by intensity redistribution, 
amplitude dependent phase-shift and relative separation distances. We denote this kind of 
collision scenario as type-II shape changing collision (type-II SCC). The most important 
consequence of type-II SCC is the possibility of soliton amplification in all the components. 
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Mixed CNLS equations are not only of mathematical interest but also of considerable 
physical significance. In particular, mixed 2-CNLS equations arise as governing equations 
? or electromagnetic pulse propagation in left handed materials with Kerr-type nonlinearity 
9| and in the modified Hubbard model in the long- wavelength approximation [llj. It can 
also be noticed that in the mixed 2-CNLS system, if we assume the fields qi and q\ (see 
Eq. (1) below) propagate in the anomalous and normal dispersion regimes, respectively, the 
self phase modulation (SPM) coefficients are positive and cross phase modulation (XPM) 
coefficients are negative in both the components. This kind of nonlinearities can be realized 
in the quadratic nonlinear materials with large phase-matching [25]. The main aim of this 
paper is to analyse type-II SCC behaviour of bright solitons in a possible integrable extension 
of mixed 2-CNLS equations involving the linear self and cross coupling terms which can be 
of physical interest. This study reveals the fact that in this system during the type-II SCC 
process there also occurs periodic energy switching due to the linear coupling terms which 
can be suppressed/enhanced by type-II SCC. The distinct feature is that the system exhibits 
large periodic energy switching between the components with very small linear self coupling 
strength. 

The plan of the paper is as follows. In Sec. II we present the statement of the problem. 
Bright soliton solutions of the mixed 2-CNLS equations with linear self and cross coupling 
terms are obtained in Sec. III. Section IV is devoted to a detailed analysis of the role of the 
coupling terms on type-II SCC. The final section is allotted for conclusion. 

II. STATEMENT OF THE PROBLEM 

Lazarides and Tsironis [9|] have obtained the following set of governing equations for elec- 
tromagnetic pulse propagation in isotropic and homogeneous nonlinear left handed materials, 

ki,z + qi,tt + 2/i (ai|gi| 2 + cr 2 \q 2 \ 2 ) qi = 0, (la) 
iq2,z + <?2,« + 2/i (cri|gi| 2 + cr 2 |g 2 | 2 ) <?2 = 0, (lb) 

by taking the effective permittivity and effective permeability to be intensity dependent and 
following a reductive perturbational approach . Here qi and q 2 are the electric and magnetic 
field components of the electromagnetic pulse, respectively, the subscripts z and t denote the 
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partial derivatives with respect to normalized distance and retarded time respectively, while 
p is the measure of nonlinearity, o\ and o"2 can be either +1 or -1. From a mathematical point 
of view the above equation reduces to the integrable mixed CNLS equation presented in Ref. 



13j when o\ = 1 and 02 = — 1. In Ref. [24] singular and nonsingular bright soliton solutions 
of Eq. (1) have been obtained. There it has been shown that even though system (1) admits 
SCC as in Manakov system, the intensity redistribution occurs in a completely different way 
which is not possible in the Manakov system [6]. As pointed out in the introduction we 
denote such a collision picture as type-II SCC. A typical type-II SCC where enhancement 
(suppression) of intensity occurs in both the components of soliton S\{S-i) after collision is 
shown in Fig. 1 (All the quantities in this and rest of the figures are dimensionless) . Note 
that the reverse scenario is also possible. 

The next natural step is to study how the linear self coupling resulting from the same 
component and the cross coupling arising from the other component influence the type- 
II SCC process. In this regard we consider the following non-dimensional mixed CNLS 
equations with linear self and cross couplings. 

iqu + Qut + pqi - Xl2 + 2// (|gi| 2 - |g 2 | 2 ) <?i = 0, (2a) 
iq2z + qut - pq2 + xqi + 2 /i (|<?i| 2 - |<?2| 2 ) q-i = 0, (2b) 

where p and x are the se lf an d cross coupling coefficients, respectively. The above system 
reduces to the Lindner-Fedyanin system [ll| with x = 0, which is a ID continuum limit of 
2D Hubbard model. We also believe that this mathematical model can be experimentally 
realized in fiber couplers fabricated from suitably engineered left handed materials. Another 
physical model associated with Eq. (2) is two component BECs coupled with two photon 
microwave field where the signs of co-efficients of SPM and XPM can be tuned suitably 
through Feshbach resonance. In fact, this type of linear couplings can be introduced in 
BECs by external microwave or radio frequency that induces Rabi or Josephson oscillation 
between population of two states [26]. Also at this point it is worth to point out that a 



similar set of equations with 
in twisted birefringent fibers 



bcusing nonlinearity (\qi\ + \q%\ ) arises in pulse propagation 
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FIG. 1: Shape changing collision of bright solitons in mixed 2-CNLS system. 

III. BRIGHT SOLITON SOLUTIONS 

To obtain the bright soliton solutions of Eq. (2), we identify a transformation which 
reduces Eq. (2) to the well known integrable mixed CNLS equations (1) with o\ — 1 and 
o"2 = —1. The transformation can be written as 



»1 









cosh(|)e irz sinh(f)e- ir2 
sinh(f)e ir2 cosh(f)e- in 





( qim | 




I q2m 1 



(3a) 



5 



where the real parameters 6 and V are expressed in terms of the coupling parameters p and 

X as 

fl = tanh- 1 ^V r =v ^^, x <p. (3b) 

In Eq. (I3al) . qi and g2 satisfy Eq. (2) with linear self and cross coupling terms included 
while qi m and q 2m satisfy Eq. (1) in the absence of linear couplings ( p = x — 0). It is 
obvious that if the cross coupling term \ becomes zero the above solution (qi, q 2 ) is the same 
as that of the mixed CNLS equations (1), with o\ = —o 2 = 1, except for a multiplicative 
phase factor e lVz in the q\ component and e~ lTz in the q 2 component. We now confine our 
analysis to the cases where (qi m , <?2m) correspond to bright soliton solutions and analyse the 
nature of {qi,q 2 ) satisfying Eq. (2) through the relation Q. 



A. Bright one-soliton solution 



With the knowledge of the bright one soliton solution of the integrable mixed 2-CNLS 
equations for (qi m , q 2m ), given in Eq. (6) of Ref. |24|, we write down the one soliton solution 
of Eq. (2) by using the transformation ([3]) as 







I*; 





cosh(|je 
sinh(|)e 



iVz 
Vz 



sinh(|)e" 
cosh(|)e 



-iVz 
-iVz 









\A 2 



/. i/,'S(xh I ViR + ^j e iriu , (4a) 



where 



A, 



h(t + ikiz) = k 1R (t - 2kuz) + i(k u t + (k\ R - k^z) = r) 1R + irjii, (4b) 

J = 1,2, (4c) 



a? 



.,Wi 



p [<Ji\a{ '\* + a 2 \a 1 



(2) I 



I (2) i 2 



1 1/2 : 



0"! 



"02 = 1- 



(4d) 



(h + k\f 

In the above expressions the suffices R and I denote the real and imaginary parts, respec- 
tively. The one soliton solution (4) for (qi,q 2 ) is characterized by three arbitrary complex 



parameters a 



(1) « (2) 

1 > "1 ; 



and k\, in addition to the real coupling parameters p and x- Also 



note that the value of \ is restricted by Eq. (3b) as |x| < \p\ since tanh^ = x/P- As in 
the case of mixed CNLS equations, solution (4) can be both singular and nonsingular. The 



condition for non-singular solution is given by \a\ 



Mi 



> a 



(2) | 



1. In this work we deal with 



nonsingular solutions only as they are of specific physical interest. 
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FIG. 2: Intensity plots of bright one soliton solution (4): (a) in the presence of coupling (p, % > 0), 
(b) in the absence of coupling (p = % = 0). 

1. Analysis on bright one-soliton solution 

Typical plot of non-singular bright one soliton solution (4) of Eq. (2) with the condition 
|x| < \p\ is shown in Fig. 2. From the figure we observe that the role of the linear coupling 
terms in (2) is to induce spatially periodic intensity switching between the two components 
qi and The periodic oscillations Fig. 2(a) during the intensity switching depends partic- 
ularly on the difference between the self and cross coupling terms (p and x) m addition to 
the soliton parameters ki, and af 1 . For comparison we have plotted the corresponding 
one-soliton solution in the absence of coupling terms in Fig. 2(b). It is interesting to note 
that this periodic intensity switching can be completely suppressed by suitably choosing A 2 

(2) 

or«j . To see this, we compute the intensity of the soliton in the two components and write 
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them as 



kf R |Ai| 2 cosh 2 



+2Ui|U 2 | cosh 



k{ R I A I sinh 2 



+ \A 2 
sinh 



sinh 



cos(2r^ + Q) 



(5a) 



+ \A 2 \ cosh 2 



+2|Ai||A 2 | cosh - sinh - cos(2r^ + Q) 



(5b) 



where 



P 
Q 



sech ( ki R (t — 2kuz) + 
An 



tan 



.4 



1H 



— tan~ 



R 
~2 

A 2 i 
A 2R 



(5c) 
(5d) 



It is clear from the above expressions that the oscillatory term cos(2rz + Q) appearing 
in (5a) and (5b) leads to the periodic oscillations during energy switching. One can also 
verify that the spatial period of oscillation is Z = ^. Thus for larger T the width of 
spatial oscillations is smaller. Also the amplitude of oscillation (21^4x1 |A 2 | cosh (|) sinh (|)) 
increases with decreasing T due to the dependence of 6 on \ an d p (see Eq. (3b)) . We also 
note from Eqs. (5a) and (5b) that the oscillatory term (third term on the right hand side) 



vanishes when \A 2 \ = 0, that is a\ = 0, or \A\\ = 0, that is a[' =0. At a first glance, it 
seems that the periodic energy switching scenario is similar to that of the Manakov system 
6] with linear coupling terms arising in the context of twisted birefringent fibers j^]. But 
the way in which the switching occurs is different due to the hyperbolic terms in Eqs. (5a) 
and (5b), since < sinh 2 8/2 < 00, < cosh 2 9/2 < 00. To be more precise, the amplitude of 
periodic oscillations and periodic switching of energy between the two components can vary 
exponentially but restricted by the condition (3b). In this connection, we would like to add 
that a quite different kind of multi-scale periodic beating of intensities without switching 
and of different physical origin as compared with the present mixed CNLS case, has been 
observed during the propagation of multisoliton complexes in the integrable N-CNLS system 
with focusing type nonlinearities 28 1. 



(i) 
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B. Two soliton solution 



The bright two soliton solution of Eq. (2) can be obtained by applying the transformation 
(3) to the two soliton solution of the integrable mixed CNLS equation given by Eq. (10) in 
Ref. [241]. The explicit form of the solution is 

q x = 1 [ f a fV r2 cosh(^) + afV ir2 sinh(^ e m 



D VA 2' 1 y 2\ 

+ (a^ 2 cosh(^) + 4 2) e- ir ^sinh(^)J e m 

+ (V 11 cosh(^)e* r2 + e Su sinh(^) e - ir *J e Vl+ ^ +m 

+ (V 21 cosh(^)e ir " + e 522 sinh(^)e- ir ^ e m+ ^ +m ^j 

q 2 = -^^«S 1) e jr2 sinh(^) + a! 2) e-^cosh(^)^e ? ' 1 

+ ^V rz sinh(~) +4 2) e-^cosh(M e m 

+ (V 11 sinh(^)e 4rz + e Sl2 cosh(^)e~ ir ^ e" 1+ ^ +r?2 



Various quantities found in Eq. (6) are defined as below following Ref. [241 ] : 



rji = ki(t + ikiz), e 6 ° = Ku - 



R, K ll R 9 K 22 

e 1 = — , e 2 — 



fa + k* fa + k{ fa + k*' 

(fa - k 2 )(a ( f ) K 21 - a { 2 j) Kii) S2 . _ {h - fa)(a^K W ~ a? k 22 ) 



. e 02 J 



(fa + fa) (fa + fa) ' (fa + k* 2 )(fa + fa 2 ) 



eRs = Ti — , , 7Jm; — ^^(^11^22-^12^21), (6d) 



\fa - fa 



|2 



(fa + kl)(k 2 + k* 2 )\fa + k 



*\2 
2 I 



and 



1 (1) (1)* 1 (2) (2)* 



(6a) 



+ [e d21 sinh(-)e ii 2 + e d22 cosh(-) e -* i 2 J , (6b) 

(i) 

where a\ s are complex parameters and the denominator D is given by 

D — 1 _|_ e m+v*+Ri _|_ e vi+V2+ s o _|_ e v!+V2+$o _|_ e V2+v2+ R 2 _|_ e m+rii+V2+ri2+ R 3 _ (g c ) 



= — -. r — , i, j = 1,2, (6e) 

3 {fa + k *j) 

where 0\ = 1 and o 2 = — 1. This solution represents the interaction of two bright solitons 
in the presence of self and cross coupling terms. Although the above solution features both 
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singular and nonsingular solutions in the following we consider only the nonsingular soliton 



solution which results for the choice 



24) 



+ 



Kn > 0, k 2 2 > 0, K11K22 - > 0, 

\ki - k 2 \ 



K\\K 22 



K 12 



2 y kmkvR 2\ki + k 



> 



(7a) 
(7b) 



"21 V n lK r »2K riTft2| 

and analyse their collision behaviour. In a similar way the multi-soliton solution of Eq. (2) 
can be obtained by applying the transformation to the multi-soliton solution given in the 
appendix of [24] with N = 2. 



IV. SHAPE CHANGING COLLISION OF SOLITONS WITH PERIODIC EN- 
ERGY SWITCHING 

We have already explained the nature of type-II SCC of bright solitons in Sec. II. In this 
section, we analyse the influence of linear cross coupling terms on the above mentioned type- 
II SCC. We perform an asymptotic analysis [24J for the choice k 1R , k 2R > and ku > k 2 j. 
To facilitate the understanding of the collision dynamics we consider the intensities of the 
two colliding solitons in the asymptotic limits at z — ► — 00 (before collision) and z —> oc 
(after collision). In their explicit forms the intensities of solitons as z — > ±00 read as 
2 

!^2 / I An: 
v nR \ A 1 



if 



P 

* n. 



kl„( |A?±| 2 cosh 2 ( " 1 +|^ ± | 2 sinh 2 



+2(^11^1 cosh - sinh - cos(2rz + Q n 



(8a) 



where 



sech (k nR (t - 2k nI z + R n )) , 



Qr 



tan 



A n± 
{hi 

A n± 
n lR 



tan 



-1 



A 



n± 
21 



An± 
A 2R 



n 



l,2,j = l,2. 



(8b) 
(8c) 



Here the quantities A n - k nR and A™ + k nR , j,n = 1,2, are the amplitudes of the nth soliton 
in jth component before and after interaction respectively in the absence of linear coupling 
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(x = P = 0)? where A n, s take the following forms before and after interaction 241 ]: 



A\~ 

A\~ 
Af 

aV 



Af 
Af 



Various quantities occuring in Eqs. 
verified that 

2 



a 



(i) 



-Rill 



(2) 



a 



(i) 
2 

(2) 



e -{R 1 +R 3 )/2 

{k 2 + k* 2 ) ' 

e -(i?2+R 3 )/2 

(ki + kl) ' 
P-R2/2 



(9a) 
(9b) 
(9c) 
(9d) 



ar I ( k 2 + k* 2 ) 
are defined in Eq. (6). From Eqs. (8) it can also be 





2 










Pn 




Pn 



k 2 nR (\A^\ 2 -\Ar\ 2 ) 



[I 



n 



1,2, 



(10) 



which is a consequence of the conservation law of intensities in the mixed CNLS system. 
The role of linear coupling parameters on type-II SCC and vice- versa can be well understood 
by analysing the asymptotic expressions (8) which clearly shows that these terms induce 
periodic switching of intensity between the two colliding solitons in both the components 
qi and q 2 - At a first sight it seems that the periodic intensity switching in a given soliton 
(say soliton Si) is influenced only by the same soliton present in the other component. But 
a careful analysis shows that the presence of the second soliton (say soliton S2) plays an 
indirect but predominant role in controlling the switching process through type-II SCC and 
vice-versa. Various possibilities of such collision scenario are given below: 

1. The coupling results in periodic oscillations in the energy switching process throughout 
the collision process due to the oscillatory term cos(2r^ + Q n ) appearing in Eq. (8a). 
As in the case of one soliton solution here also the amplitude and width of the periodic 
oscillations increase with decreasing V. Thus the important feature of such collision 
process is that the amplitude of periodic energy switching can be large depending upon 
the relative signs of linear coupling terms p and x- This periodic energy switching 
behaviour, in the presence of coupling, depends on the a and fc-parameters and also 
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FIG. 3: Type-II shape changing collisions with periodic intensity switching. 



on the linear coupling coefficients. Thus the oscillating energy switching process co- 

a (i) Q (2) 

exists with type-II SCC for -4^ ^ -fey. Such a two soliton collision process with 
periodic intensity switching is shown in Fig. 3 for = 0.7226 + 1.12542, ot^ = 
0.8484 + 0.2625i, a£ ] = 0.5511 + 0.8584i, a { p = 0.1923 + 0.0595i, p = 1, x = 0.5, 
ki = 1 + k 2 = 1.1 — i. Eq. (8) also shows that the coupling enhances the amplitude of 
the soliton in a given component before and after interaction due to the contribution 
from the other component as compared with the bright soliton collision case in the 
absence of coupling. 

2. The distinct feature of this collision process is that the intensity redistribution can 
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FIG. 4: Suppression of periodic oscillations in 52 after interaction in a type-II SCC process. 

be used to control the switching dynamics. One interesting possibility is complete 
suppression of oscillation either before or after collision in a particular soliton say 
u S n " by making any one of \A™~\ or j, n = 1,2, to be zero, respectively, 

with commensurate changes in the other soliton. As the nonsingular condition (7a) 
of the solution rules out the possibility of making \A™ ± \ to be zero , the complete 
suppression of periodic oscillation of intensities in both the components of soliton 
S n before (after) collision can be obtained by choosing \A%~\ (\A2 + \) = 0. This 
suppression (enhancement) of intensities of a particular soliton in a given component 
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FIG. 5: Suppression of periodic oscillations in S\ before interaction in a type-II SCC process. 

during the type-II SCC results in the enhancement (suppression) of amplitude of 
periodic oscillations in the other colliding soliton as inferred from Eq. (8). Fig. 4. 
shows the type-II SCC scenario in which the oscillations in the q\ and q 2 components 
of S 2 after interaction are completely suppressed, for the choice 6093+0. 9489i, 

= 0.4978 + 0.1540i, 4 1) =0.5403+0.8415i, a { 2 2) = 0, p = 1, x = 0.5, h = 1 + i, 
k 2 — 1.1 — %■ The reason for this is that in the absence of coupling terms soliton S 2 
undergoes type-II SCC with Si and its intensity in q 2 component after interaction is 
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FIG. 6: Elastic collision of bright solitons with periodic oscillations in mixed CNLS system with 
linear self and cross couplings. 



exactly zero for the given parametric choice. Similarly Fig. 5 shows the suppression of 
oscillations in q\ and q 2 components of S\ before interaction for the parametric choice 
= 1, af ] = 0, c4 1} = 1.0201, a { 2 2) = 0.2013, p = 1, x = 0.5, h = l + i, k 2 = 1.1 -i. 
This kind of switching process arises from the fact that in the absence of coupling the 
intensity of S\ in q 2 component (that is, \A\~\ 2 k\ R ) is zero before it collides with S 2 - 

3. The standard elastic collision with periodic energy switching only arises for the choice 

a (l) a (2) (1) 

-4jt = -fey. This is shown in Fig. 6 for the parametric choice a\ ' = 0.6782 + 1.0562i, 
af = 0.6782 + 1.0562i, = 0.7247 + 0.2242i, a { 2 2) = 0.7247 + 0.2242i, p = 1, 
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X = 0.5, ki—1 + i, k 2 = 1.1 — i. 

In order to appreciate the significance of the present system, we compare the soliton collision 
behaviour with that of twisted birefringent fibers 2jj] which involve focusing type nonlin- 
earities. The crucial difference follows from Eq. (9), which says that the energy exchange 
between the two components (51,52) before and after collision is constant and as a result a 
given soliton experiences the same effect (either suppression or enhancement of intensity) 
in both the components during its collision with other soliton contrary to the twisted bire- 
fringent system. Thus the amplitude of oscillation due to coupling can be simultaneously 
enhanced/suppressed after collision in both the components as a consequence of type II- 
SCC, a situation which is not possible in twisted birefringent fibers. Another important 
advantage is the efficiency of switching due to linear couplings. Here the coupling terms 
influence the energy switching exponentially due to the hyperbolic terms (see Eq. (9)). 
This suggests large switching of energy with small self coupling strengths, as compared with 
Manakov system with linear couplings, a desirable property in fiber couplers. 



V. CONCLUSION 



In this paper, we have shown that the set of mixed 2-CNLS equations with linear self 
and cross coupling terms can be transformed to the standard integrable mixed 2-CNLS 
equations by performing the transformation (3a). The bright soliton solutions are obtained 
by applying this transformation to the recently reported bright soliton solutions of the mixed 
2-CNLS equations [24j without linear coupling terms. Our study shows that inclusion of 
linear self and cross coupling terms lead to periodic energy switching among the components. 
We have also pointed out that due to the exponential dependence on the coupling terms, 
the energy switching can be large with small coupling strengths. In a two soliton collision 
process such periodic energy switching coexists with the type-II shape changing collision 
behaviour. However the standard elastic collision process can take place with or without 
periodic energy switching for very specific parametric choices. An important result which 
follows from the present study is that the shape changing collision of type-II can be used 
suitably to suppress or enhance the periodic oscillations in the energy switching process 
completely or partially and also simultaneously in both the components. These results can 
give further impetus in understanding the Lindner- Fedyanin system in the continuum limit, 
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and can also find potential applications in fiber couplers and in BECs. 
ACKNOWLEDGMENTS 

T. K acknowledges the support of Department of Science and Technology, Government 
of India under the DST Fast Track Project for young scientists. He is also grateful to the 
constant support of the Principal and Management of Bishop Heber College, Tiruchirapalli. 
The work of M. V and M. L are supported by a DST-IRPHA project. M. L is also supported 
by a DAE-BRNS Raja Ramanna Fellowship. 



[1] N. J. Zabusky and M. D. Kruskal, Phys. Rev. Lett. 15, 240 (1965). 

[2] Yuri S. Kivshar and Govind P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals, 

(Academic Press, San Diego, 2003). 
[3] M. Lakshmanan and S. Rajasekar, Nonlinear Dynamics: Integrability, Chaos and Patterns, 

( Springer- Verlag, Berlin, 2003). 
[4] N. Akhmediev and A. Ankiewicz, Solitons: Nonlinear Pulses and Beams, (Chapman and Hall, 

London, 1997). 
[5] A. C. Scott, Phys. Scr. 29, 279 (1984). 

[6] S. V. Manakov, Zh. Eksp. Teor. Fiz. 65, 505 (1973) [Sov. Phys. JETP 38, 248 (1974)]. 
[7] C. R. Menyuk, IEEE J. Quantum Electron. 23, 174 (1989). 

[8] D. N. Christodoulides, T. H. Coskun, M. Mitchell, and M. Segev, Phys. Rev. Lett. 78, 646 
(1997). 

[9] N. Lazarides and G. P. Tsironis, Phys. Rev. E 71, 036614 (2005). 
[10] B. P. Anderson and M. A. Kasevich, Science 282, 1686 (1998). 

[11] V. G. Makhankov, Phys. Lett. A 81, 156 (1981); U. Lindner and V. Fedyanin, Phys. Stat. 

Sol. (b) 89, 123 (1978). 
[12] V. E. Zakharov and E. I. Schulman, Physica D 4, 270 (1982). 

[13] V. G. Makhankov, N. V. Makhaldiani, and O. K. Pashaev, Phys. Lett. 81A, 161 (1981) ; V. 
G. Makhankov and O. K. Pashaev, Teor. Mat. Fiz. 53, 55 (1982). 



17 



[14] R. Radhakrishnan, R. Sahadevan, and M. Lakshmanan, Chaos, Solitons and Fractals 5, 2315 
(1995). 



[15 
[16 

[17; 

[18 
[19 
[20 
[21 
[22 
[23 
[24 

[25 
[26 

[27; 

[28 



R. Radhakrishnan, M. Lakshmanan, and J. Hietarinta, Phys. Rev. E 56, 2213 (1997). 
T. Kanna and M. Lakshmanan, Phys. Rev. Lett. 86, 5043 (2001) . 
T. Kanna and M. Lakshmanan, Phys. Rev. E 67 , 046617 (2003). 
Darren Rand et al, Phys. Rev. Lett. 98, 053902, (2007). 
C. Anastassiou et al, Phys. Rev. Lett. 83, 2332 (1999). 

M. H. Jakubowski, K. Steiglitz, and R. Squier, Phys. Rev. E 58, 6752 (1998). 
K. Steiglitz, Phys. Rev. E 63, 016608 (2000). 

M. Lakshmanan and T. Kanna, Pramana - J. of Phys. 57, 885 (2001). 

T. Kanna, E. N. Tsoy, and N. N. Akhmediev, Phys. Lett. A 330, 224 (2004). 

T. Kanna, M. Lakshmanan, P. Tchofo Dinda, and Nail Akhmediev, Phys. Rev. E 73, 026604 

(2006). 

A. G. Kalocsai and J. W. Haus, Phys. Rev. E 52, 3166 (1995). 



B. Xiang, [arXiv:nlin.P S/0703001vl/ Mar. 2007. 

R. Radhakrishnan and M. Lakshmanan, Phys. Rev. E 60, 2317 (1999). 
Nail Akhmediev and Adrian Ankiewicz, Chaos 10, 600 (2000). 



18 



